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Abstract. The 2+1 dimensional pure SU (N) gauge theories with N < 4 are candidates for applying the powerful tools of 
scaling and universality to their deconfinement transitions at finite temperature. The corresponding 2 dimensional g-state 
Potts models with q < 4 have 2 nd order transitions, and we can exploit many exact results to obtain accurate critical couplings, 
transition temperatures, critical exponents, and the leading behavior of the continuum string tension near the phase transition 
on one side, together with its dual on the other. Thereby, the self-duality of the 2d spin models is reflected in a duality between 
spacelike vortices and confining electric fluxes. We also discuss the relevance of center symmetry and the corresponding 
vortices for confinement in full QCD when the electromagnetic interactions of fractionally charged quarks are included. 
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Introduction 

The finite temperature deconfinement transition in 
SU(N) gauge theories in d+ 1 dimensions is very well 
understood in terms of the spontaneous breakdown of 
their global Z N center symmetry [1]. Static fundamental 
charges represented by Polyakov loops P(x) transform 
under this symmetry like spins s, in a d dimensional q- 
state Potts model with q — N and Hamiltonian [2], 

J? = -J £ 8 Si>Sj - £ 5,,o , s t = 0, 1, . . .q - 1 , (1) 

with nearest neighbor coupling J. A non-zero external 
field H, inversely related to the quark mass m q , may be 
included to mimic the leading effect of heavy dynamical 
quarks. When l/m q — 0, the Polyakov loop develops 
a non-zero expectation value only in the deconfined, 
Zjv-broken phase, while the expectation value of P(x) 
vanishes in the disordered, confined phase much like the 
spontaneous magnetization in the spin model. 

This is well described in terms of spacelike center vor- 
tices which play the role of spin interfaces. They sepa- 
rate regions where the Polyakov loop differs by a phase 
z € Z N , so that their proliferation disorders the Polyakov 
loop and leads to confinement. The suppression of spa- 
tial center vortices at high temperatures coincides with 
the ordering of the Polyakov loop, and their free energy 
offers an elegant order parameter for the transition [3]. 

When the phase transition is of second order its de- 
scription is universal [4]. In 3+1 dimensions this only ap- 
plies to 5(7(2), where spatial center vortex sheets show 
the universal behavior of interfaces in the 3d Ising model 
[5]. In 2+1 dimensions, on the other hand, both SU(2) 
and SU (3) exhibit a second-order deconfinement transi- 
tion, and many exact results from the interfaces in the 2d 
Ising and 3-state Potts models can be fruitfully exploited 
for precision studies of the (2+l)d gauge theory [6, 7, 8]. 



Even for 5(7(4), which was previously found to have a 
weak first-order transition [9], at least approximate Potts 
scaling can be observed in a wide range of intermediate 
length scales near criticality [10, 8]. 

Electric Fluxes and Self-Duality in 2+ld 

Interfaces in the spin models are typically introduced 
as frustrations along which the coupling of adjacent spins 
favors cyclically shifted spin states rather than parallel 
ones for the usual ferromagnetic couplings J > 0. They 
form d — 1 dimensional surfaces dual to links at which 
the 8 SUS j in Eq. (1) are replaced by 8 ShS . +m mo d q , and are 
conveniently studied by introducing analogous cyclically 
shifted boundary conditions. 

The ratios R q m ' n) = Z q m ' n) /zf fi) of g-state Potts model 
partition functions Z q with combinations of cyclically 
shifted (ra,n)-boundary conditions on a 2d torus over the 
periodic ensemble with m = n = then define the inter- 
face free energies per temperature as Fj m,nS> — — lnR q m,n ^ . 

Center vortex free energies in 7+1 dimensional SU(N) 
are defined from analogous ratios of partition functions 
with 't Hooft's twisted boundary conditions [11] over the 
periodic ensemble. These are classified either as mag- 
netic twists, defined in purely spatial planes, or as tem- 
poral twists in the planes oriented along the Euclidean 
time direction. The former are irrelevant for the phase 
transition [12] and will not be further discussed here. The 
latter are labeled by kj ■■= 0, 1 , . . . N— 1 , and introduce spa- 
tial center vortices perpendicular to the spatial direction 
;' of the twist. These separate regions where fundamental 
Polyakov loops differ by a center phase z — e 2mk '/ N and 
are thus like the spin interfaces in the Potts models. 

Exact maps between the spin systems and their dual 
theories, in terms of disorder variables on the dual lattice, 
are provided by Kramers-Wannier duality [13]. For the 
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FIGURE 1. The ratios of St/ (2) partition functions R k (x) 
for twist (1,0) and R e (x) for electric flux (1,0) over the FSS 
variable x oc ±L/| (data here for N t =4 and up to 96). 

g-state Potts models in 2 dimensions, these dual theories 
are g-state Potts models again, but at a dual temperature 
T which is swapped around criticality at T c . A simple 
proof for infinite lattices is given in [2]. Duality transfor- 
mations in a finite volume do not preserve boundary con- 
ditions, however. Periodic boundary conditions on one 
side generally correspond to fluctuating boundary condi- 
tions on the other. From the exact finite-volume duality 
transformation for the 2d g-state Potts models [7], 

y f{rn-sm) 

\k) = Us 



L, S <"V) 

where the coupling per temperature K = J/T and its 
dual K = J/f are related by (e* -l)(e K -l) = q, with 
criticality at K = K = K c = ln(l + y/q). 

With q = N, the discrete 2d Z^-Fourier transform (2) 
precisely resembles 't Hooft's relation between the ratios 
Rk(k) = Zk(k)/Zk(0), with temporally twisted k b.c.'s 
over the periodic ensemble, and those of electric fluxes e 
relative to the no-flux ensemble with fluctuating temporal 
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In the (2+l)d gauge theory, temperature is the same on 
both sides of the Z^-Fourier transform. As a consequence 
of the self-duality of the spin model, however, the free 
energies of spatial center vortices and those of the con- 
fining electric fluxes are mirror images of one another 
within the universal scaling window around a second or- 
der phase transition. This is shown for SU (2) in Fig. 1, 
where rolin the spin model amounts to x o —x, and 
the same is true for 51/(3) [8]. 

Critical Couplings and Finite-Size Scaling 

The by far most efficient method to determine critical 
couplings j5 c from simulations in finite volumes is based 
on self-duality [8]: with that one has R e (P) = Rk(P) for 




FIGURE 2. Critical couplings for (2+l)d St/ (2) (N t = 4) 
from self-duality compared to those of Ref. [6], with N s 
up to 96 and the infinite-volume extrapolated result j8 c = 
6.53661(13), shown as the narrow grey band here. 

matching e and k at lattice coupling j5 = j3 c . Defining 
pseudo-critical couplings in a finite volume by this re- 
quirement, R e (P) =Rk(P)> it is straightforward to show 
that the leading finite-size corrections do not change its 
value. This is because they give the same j3 -independent 
contribution to both, R e and R^, and hence move their 
intersection point vertically, without any shift in j3. 

In Fig. 2 we demonstrate how well this works by 
comparing the results from self-duality in (2+l)d SU(2) 
to the high-precision determination of [6]. The method 
used there allowed a finite-size scaling ansatz to fit the 
pseudo-critical couplings shown as the lower points in 
Fig. 2. Their infinite-volume extrapolation led to the 
value /3 £ = 6.53661(13) for N t = 4 time slices and spa- 
tial lattice sizes up to N s — 96, which is indicated to- 
gether with its small error by the grey band in the figure. 
Within the present statistical errors, the pseudo-critical 
couplings from self-duality are fully consistent with this 
narrow band from the previous (far more expensive) de- 
termination already for N s = 16, i.e., we practically re- 
produce the infinite volume result with an aspect ratio of 
only 4:1. The weighted mean from N s = 16,24 and 32 
yields J3 C = 6.53651(13). 

With accurate critical couplings we can assess the 
finite-size scaling near the second order phase transition 
at t = T /T c — 1 = 0, where generalized couplings such 
as the vortex-ensemble ratios R^, for sufficiently large 
sizes L, only depend on L 1//y f in a universal way. This 
can be used, for example, to extract the correlation- 
length critical exponent v from the slope s of the center 
vortex free energies /^(/3) = — InR^ in j3 at j8 c . Because 
t oe (J3 — j3 c ), the slope grows with the spatial lattice size 

as s ~ N^ V . Our present data for (2+l)d SU(N) in this 
way yields v = 0.99(2), 0.82(4) and 0.60(5) for N = 2, 3 
and 4, as compared to v = 1,5/6 and 2/3 for the 2d Potts 
models with g = 2, 3 and 4, respectively. 

Moreover, for St/ (2) this dependence of the center 
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FIGURE 3. SU (2) center- vortex partition functions with A',- 
dependent rescaling of x for N, = 4 to 10 and N s up to 96. 

vortex free energies is known exactly, from the analyt- 
ically computable universal scaling functions Fj(x) for 
the interface free energies in the 2d Ising model [7, 14]. 
One obtains Fk(x) = Fj(-Xx) where the single non- 
universal factor A can be accurately determined from 
one-parameter fits near x — LT c t = 0. This was done for 
N t = 4 to 10 in [7], which allows the continuum extrap- 
olation A = 1.354(25) for Nt — > °°. We can furthermore 
compensate discretization effects by simply rescaling the 
finite-size scaling variables on different N t lattices as 
x — s- X(N t )x, to observe good (continuum) scaling around 
criticality for all available N s and N t , as shown in Fig. 3. 

Other immediate consequences of this extrapolation 
are, e.g., the behavior of the continuum string tension and 
its dual around the deconfinement transition in (2+l)d 
SU{2): o/T c 2 = 2.387(44) \t\ + ■■■, for t -> 0", and 
a/T c = 2.387(44) t + ■ ■ ■ , for t -> 0+, respectively [7]. 

Conclusions and Outlook 

Universality and scaling allow one to study in de- 
tail the deconfinement transition in the 2+ 1 dimensional 
SU(2), 51/(3), and at least approximately also SU(4) 
gauge theories, by relating their spatial center vortices 
to interfaces in 2d Potts models for which many exact 
results are available. The self-duality of these models is 
reflected in the gauge theories: around criticality, the free 
energies of the confining electric fluxes are mirror im- 
ages of those of spatial center vortices. We demonstrated 
how this can be exploited to remove the leading finite- 
size corrections in the determination of critical couplings 
from numerical simulations. These then allow detailed 
finite-size scaling analyses and simple continuum extrap- 
olations such that the consequences of universal scaling 
are carried over to the continuum gauge theory. 

We close with a short remark on the relevance of cen- 
ter symmetry and center vortices for confinement in full 
QCD when including the electromagnetic interactions of 
fractionally charged quarks. Dynamical quarks explicitly 
break center symmetry in the same way as the external 
field H in (1). When electromagnetic C/(l) couplings are 
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FIGURE 4. 51/ (2) Polyakov loop with half-integer charged 
dynamical quarks ( K = 0. 1 5) on an 8 3 x 4 lattice with hot 1/(1) 
start to provide Z2 disorder in Coulomb phase (at large j8 em ). 

included, however, this may become more like a fluc- 
tuating external field. Because of the quarks' fractional 
charges, the 1/(1) gauge action is then unable to remove 
initial disorder as it evolves from strong to weak coupling 
for integer electric charges. This is demonstrated for a 
ST/ (2) x U(\)/Z2 toy model with half-integer charged 
quarks in Fig. 4, where the quenched behavior is recov- 
ered for the SU (2) Polyakov loop with a hot U(l) start, 
far into the Coulomb phase for integer charges, see [15]. 
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